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;-H ■ ABSTRACT 

, Aims. We analyse the point- symmetric planetary nebula NGC 6309 in terms of its three-dimensional structure and of internal varia- 

tions of the physical conditions to deduce the physical processes involved in its formation. 

Methods. We used VLA-D /l3.6-cm continuum, ground-based, and HST- archive imaging as well as long slit high- and low-dispersion 
spectroscopy. 

Results. The low-dispersion spectra indicate a high excitation nebula, with low to medium variations of its internal physical conditions 
. (10,600K < re[Oiii] < 10,900 K; 10, lOOK < re[Nii] < 11, 800 K; 1440cm-3 < A^e[Sii] < 4000cm-^ HOOcm'^ < A^e[Cliii]< 

■ 2600 cm~^; 1000 cm~^ < A^e[Ariv] < 1700 cm~^). The radio continuum emission indicates a mean electron density of ^ 1900 cm~^, 
' emission measure of 5.1 x 10^ pc cm~^, and an ionised mass M(Hii)^ 0.07 M©. In the optical images, the point- symmetric knots show 
I a lack of [N ii] emission as compared with similar features previously known in other PNe. A rich internal structure of the central 

region is seen in the HST images, resembling a deformed torus. Long slit high-dispersion spectra reveal a complex kinematics in 

■ the central region, with internal expansion velocities ranging from ^ 20 to 30kms~^ In addition, the spectral line profiles from the 
QQ ' external regions of NGC 6309 indicate expanding lobes 40kms~0 as those generally found in bipolar nebulae. Finally, we have 

I found evidence for the presence of a faint halo, possibly related to the envelope of the AGB-star progenitor. 

Conclusions. Our data indicate that NGC 6309 is a quadrupolar nebula with two pairs of bipolar lobes whose axes are oriented PA=40° 
and PA=76°. Equatorial and polar velocities for these two pairs of lobes are 29 and 86 km s~^ for the bipolar system at PA=40° and 
" 25 and 75 km for the bipolar system at PA=76°. There is also a central torus that is expanding at 25 km s~^ Kinematical age for 

rS ! all these structures is around 3700 to 4000 yr. We conclude that NGC 6309 was formed by a set of well-collimated bipolar outflows 

(jets), which were ejected in the initial stages of its formation as a planetary nebula. These jets carved the bipolar lobes in the previous 
AGB wind and their remnants are now observed as the point- symmetric knots tracing the edges of the lobes. 

Key words, planetary nebulae: individual: NGC 6309 - ISM: kinematics - ISM: abundances 



1. Introduction 

In spite of all the work that up to now has been done to under- 
stand the origin of planetary nebulae (PNe) morphologies and 
their evolution (e.g., Kwok, Purton & FitzGerald 1978; Kahn & 
Westtt985i Balick^^ Icke [19881 M ellema ll995i Perinotto 
et al. 12004: Rijkhorst, Mell ema & Icke l2QQ5l Schonberner et 
al. 120051 Schonberner et al. I20Q7K some morphological struc- 
tures still remain as unsolved problems. In particular, the so- 
called point- symmetry (Stanghellini, Corradi & Schwarz, 1993; 
Gon9alves et al. 2003 ) appears enigmatic. Some observational 
studies about point- symmetric PNe (e.g. Miranda & Solf 1992> 
Lopez, Meaburn & Palmer [T993I ) relate the formation of these 
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objects to collimated outflows from a precessing central source 
known as bipolar rotating episodic jets (BRETs, see Lopez, 
Vazquez & Rodriguez 1995 ). However, in some cases, there are 
not definite proofs about the jet nature of these features as they 
are not confirmed by spectroscopic studies (e.g. Vazquez et al. 
I1999al Vazquez et al j2QQ2l) . Although some theoretical models 
have intended to explain point- symmetric PNe (e.g., Cliff'e et 
al.[T995t Livio & Pringle 1996, 19^97; Garcia-Segura & Lopez 
2000; Rijkhorst, Icke & Mellema 120041) , the origin and shaping 
of this kind of objects remains a puzzle. 

NGC 6309 is a PN whose morphology strongly suggests a 
BRET scenario for its origin, as can be seen in the HQf+[N ii] and 
[O III] images by Schwarz, Corradi & Melnick (11992b. It has two 
prominent, point- symmetric 'arms' formed by pairs of conden- 
sations in addition to a bright internal elliptical structure. Some 
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Fig. 1. Images of NGC 6309 in the light of Ha (top left), [Nii] (top right) and [Oiii] (bottom left), including an unsharp masking of 
the [Oiii] image (bottom right). North is up and east is left in all the panels. Small frames are included inside the upper panels to 
show the structure of the central region. 



previous studies on NGC 6309 have determined an expansion 
velocity Vexp[Oiii]=34kms"^ (Sabbadin I 19841 ), as well as mean 
physical parameters and total abundances (Gorny et al. 12004), 
namely, electron density and temperature A^e[S ii]=2600cm"^, 
Te [N II] =1 2,097 K, re[0 III] =1 1,845 K; and abundance ratios 
He/H=0.1, N/H=8.20, 0/H=8.64, Ne/H=7.82, and S/H=6.49 
(log H= 12, except for He/H ratio). Armour & Kingsburgh ( 2001b 
found similar values. NGC 6309 has also been included in the list 
of PNe with low ionisation structures (LIS) by Gonial ves et al. 
(120031) . Finally, the central star of NGC 6309 has been classified 
as a "weak emission line star" by Gorny et al. ([20041) . 

In spite of this suggestive morphology, an internal kinematic 
study of NGC 6309, as well as a full analysis of the physical 
conditions in the diff'erent nebular regions, has not yet been 



done. In this paper, we have carried out a radio-optical study of 
NGC 6309, including radio continuum mapping, optical imag- 
ing, and long-slit optical spectroscopy, in both, high- and low- 
dispersion. We discuss our results on the morphology, kinemat- 
ics, physical conditions, ionic abundances, and nature of the gas 
emission to explore the possible mechanisms involved in the for- 
mation of NGC 6309. 



2. Observations and results 

2.1. CCD ground-based imaging 

We obtained narrow-band direct images with the Calar Alto 
Faint Object Spectrograph (CAFOS) in the 2.2 m telescope at 
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Calar Alto Observatory (CAHAQ) in 1998 July 8. We used two 
filters centered in Ua (AA = 15A) and [Nii]6583 (AA = 20A) 
as well as a CCD Loral with 2048 x 2048 pixels. The exposure 
time for both images was 1800 s. The scale was 0^/33 pixel"^ and 
the seeing was about Vfl. We obtained an additional image on 
2004 August 3 with the 1.5 m telescope at the Observatorio de 
Sierra Nevada 0. The detector was a RoperScientific VersArray 
CCD with 2048x2048 pixels, each of a.'232 pixer^ . We used an 
[O iii]5007 (A.1 = 50A) filter with an exposure time of 900 s. The 
images were reduced following standard procedures within the 
MIDAS package. In this case, seeing was IV 3. 

Figure [T] shows the images of NGC 6309 in Ha, [Nii] and 
[O III] . An additional unsharp masking [O iii] image is also pre- 
sented in Fig. [T] to show the link between internal and external 
nebular structures. In these images, the high contrast between the 
brightness of the central region and that of the 'arms' is evident. 
In particular, the lack of emission of [Nii] from the 'arms' is un- 
expectedly remarkable given that, in general, this kind of point- 
symmetric microstructures is related to LIS (Gon9alves et al. 
1200 3 ). However, there are also some cases in which the emission 
of [Nil] is marginal (e.g. IC5217; Miranda et al. l2QQ6l) . In the 
case of NGC 6309, the lack of [Nii] from the point- symmetric 
knots indicates that these cannot be considered as LIS. 

The structure of the 'arms' consists of at least four pairs of 
point- symmetric knots, some of them clearly extended perpen- 
dicular to their corresponding radial vector from the central star. 
The central region of the NGC 6309 is a bright ellipse (major 
axis ^ 20'', PA -14°) embedded in the main body of the neb- 
ula. The arms appear to leave from the vertexes of this ellipse. 
The unsharp masking [O iii] image also shows a conelike struc- 
ture with its base located on the central ellipse and the vertex on 
the NE arm. In [O iii], faint emission connects the arms with the 
bright ellipse and seems to trace asymmetrical lobes. In fact, the 
appearance of NGC 6309 in [O iii] resembles that of a quadrupo- 
lar PN in which one of the outflows (at PA 40°) protrude into 
the lobes of the other (at PA 76°). Part of the SW lobe appears 
to be open. Finally, a circular faint halo (probably spherical) is 
detected in [O iii] image (Figs. [T] and [2]). Its center coincides with 
the central star of the nebula and its size is ^ 56''. 

In Fig. [21 we show the unsharp masking [O iii] image and 
label some morphological features that will be discussed later. 
Knots in the NE arm are named El, E2, E3, and E4, whereas 
those in the SW arm are named Wl, W2, W3, and W4. Regions 
in the ellipse are called Rl and R2 for the tips of its major axis, 
and R3 and R4 for those along the minor axis. 

2.2. HST-WFPC2 imaging 

An Hubble Space Telescope (HST) broadband image from the 
MAST Archival was used in order to improve the general view 
of the internal morphology of NGC 6309 (proposal ID: 6119; 
PI: H. E. Bond; Date of observation: 1995 August 26). Figure[3] 

^ The Centro Astronomico Hispano-Aleman (CAHA) at Calar Alto, 
is operated jointly by the Max-Planck Institut fiir Astronomie and the 
Instituto de Astrofisica de Andalucia (CSIC). 

^ The Observatorio de Sierra Nevada is operated by the Consejo 
Superior de Investigaciones Cientificas through the Instituto de 
Astrofisica de Andalucia (Granada, Spain) 

^ Some of the data presented in this paper were obtained from 
the Multimission Archive at the Space Telescope Science Institute 
(MAST). STScI is operated by the Association of Universities for 
Research in Astronomy, Inc., under NASA contract NAS5-26555. 
Support for MAST for non-HST data is provided by the NASA Office of 
Space Science via grant NAG5-7584 and by other grants and contracts. 




Fig. 3. HST broad-band CCD image of NGC 6309 in the filter 



shows the 140 sec image using the filter F555W (nearly Johnson 
V). This image can be compared with the optical ground-based 
images in Fig. [T] The bright ellipse is formed by many clumps 
as well as diff'use gas. The ellipse is open in its SW and SE re- 
gions, whereas an apparently double structure is observed in the 
S and N regions. In addition, a system of faint bubbles is detected 
towards the NE, which is probably related to the cone-like struc- 
ture seen in [O iii] (see above). Given the short exposure time of 
the HST image, evidence of the "arms" is very marginal. 

2.3. A3.6-cm VLA-D 

We obtained radio continuum observations at /l3.6-cm, toward 
NGC 6309 with the Very Large Array (VLA) of the National 
Radio Astronomy Observatory (NRAOjj in the D configura- 
tion during 1996 August 8. The standard VLA continuum mode 
with a 100 MHz bandwidth and two circular polarizations was 
employed. The flux and phase calibrators were 1331 -h 305 
(adopted flux density 5.2 Jy) and 1733 - 130 (observed flux den- 
sity 10.7 Jy), respectively. We set phase center at a(2000) = 
17^4^03^.6, ^(2000) = -12°54'37'^ The on-target integration 
time was 28 min. We calibrated the data and processed it us- 
ing standard procedures of the Astronomical Image Processing 
System (AIPS) package of the NRAO. We obtained a cleaned 
map of NGC 6309 using the ta sk IMAGR of AIPS (parame- 
ter ROBUST= -3; Briggs [T995]) . Self-calibration was also per- 
formed resulting in a final synthesized beam of 13V 9 in diameter 
and rms noise of the map of cr = 33 yuJy beam"^ 

The detected radio continuum emission is related to the 
central region and appears partially resolved. Based on a 2-D 
Gaussian fit (task IMFIT in AIPS), we estimate a size of ^ 19'' x 
16'' at PA= -27° (see Fig.©. The position of the intensity peak 
was found at Qf(2000) = 17^4^04^.28, ^(2000) = -12°54'37". 



^ The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under cooperative agreement by 
Associated Universities, Inc. 
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Fig. 2. Unsharp-masking image of NGC 6309 in the light of [O iii]. North is up and east is left. Main morphological features as well 
as some regions of particular interest are pointed out and labeled. 



The total flux density from the radio map is ^ 115-mJy, which 
is similar to those measured by Milne & Aller (11982 1): Ratag & 
Pottasch (1991); and Condon & Kaplan (19W); at 2 cm (146- 
mJy); 6 cm (102-mJy); and 21cm (132-mJy); respectively. We 
did not detect emission from the 'arms'. 



We also derived the mean physical conditions, based on 
the formalism of Mezger & Henderson (1967). Optically-thin 
emission and Tq = 10, 000 K were assumed, obtaining the 
following results: A^, = 1900 cm'^ M(Hii)=0.07 M©; EM = 
5.1 X lO^pccm"^. We assumed distance to the nebula to be 
2 kpc as the average value from the difl'erent estimates, which 
range between 1.1 and 2.5 kpc (e.g.. Daub [19821 Phillips & 
Pottasch HSU Amnuel, et al l 19541 Maciel[T98l Cahn, Kaler 
& Stanghellini,[l992|. 



2.4. Long slit low-dispersion spectroscopy 

Low-dispersion optical spectra were obtained with the Boiler 
& Chivens spectrometer in the 2.1m UNAM telescope at the 
San Pedro Martir Observatory (OAN-UNAM) in 1999 July 19- 
20 (grating of 300 lines/mm) and 2002 August 7-8 (grating of 
400 lines/mm). In afl of these cases, a CCD Tek 1024 x 1024 
was used as detector. We set the slit width to 220 yum (1.6'')- 
The spatial scale is IV 05 pixel" ^ whereas the spectral scale is 3 
and 4 A pixel"^ for the gratings of 400 and 300 lines/mm, respec- 
tively. We set the slit at several position angles, which are labeled 
as C, J, K, L (grating of 300 lines/mm), and A, C, and E (grating 
of 400 lines/mm), in Fig. [51 

Table [T] shows the dereddened spectra (considering case B of 
recombination) and physical conditions obtained with the Five- 
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Fig. 4. Contour plot /13.6-cm continuum map of NGC 6309, 
overimposed on a Ha image. Contour levels are 10, 20, 30, 40, 
50, 60, 70, 80, and 90% of the peak flux (6.3 x 10"^ Jy). The half- 
power beam width (13V 9 in diameter) is shown in the bottom-left 
corner. 




Fig. 5. Slits used in spectroscopy overimposed on [O iii] unsharp 
masking image. 

Level Atom Diagnostics Package (nebular) from IRAlQ (apply- 
ing the extinction law by Howarth 1983 for dereddening) for the 
regions Rl, R2, R3, R4, E2, Wl, and Wl. 

Physical conditions can only be determined for the re- 
gions of the central ellipse. Electron temperature and densi- 
ties obtained from the high-excitation emission lines appears 



^ The Image Reduction and Analysis Facility (IRAF) is distributed by 
the National Optical Astronomy Observatories, which are operated by 
the Association of Universities for Research in Astronomy, Inc., under 
cooperative agreement with the National Science Foundation. 



to be quite uniform (re[Oiii] = 10750K, A^e[Cliii]=2350cm-^ 
A^e[Ariv] = 1430cm"^), with variations within the estimated un- 
certainties. However, those physical parameters obtained from 
low-excitation emission ions appear as separated in two diff'erent 
regimes (re[Nii] = 10200K, and 11 800 K; A^e[S ii]=1620cm-3 
and 3800 cm"^). 

The dereddened absolute HjS flux values from regions Rl, 
R2, R3, and R4 are similar (see Table [T]). However, the cor- 
responding logarithmic extinction coefficients c^g spread from 
0.70 to 0.97 (mean value including all regions is 0.87), probably 
due to differences in the internal dust distribution, as has been 
found in other PNe (e.g. Vazquez et al. 11999b[ Lee & Kwok 
[20051) . 

In addition, we present the more common ratios for plasma 
diagnostics at the bottom of Table [T] which show that these re- 
gions are emitting by photoionisation-recombination processes 
with minimal or nonexistent shock-cooling contributions. This 
could be expected in the surroundings of a hot luminous star, 
however there are several cases of PNe in which shock-cooling 
processes make an important contribution to the microstructures 
emission (see Gonial ves 2003 for a compilation). For compar- 
ison, previous spectroscopic studies of NGC 6309 by Gorny et 
al. ('2007) give values of c(HyS) ^ 0.88, r,([Oiii])^ll 845 K, 
r,([Nii])^12097 K, and A^,([Sii])^ 2600 cm'^ 

We have also estimated ionic and elemental abundances for 
some regions using IRAF and the ICF (lonisation Correction 
Factors) method (Kingsburgh & Barlow 1994), although we are 
aware that recently, Gon9alves et al. (2006) have noticed that 
elemental abundances derived by using this method are overes- 
timated if they are obtained from narrow, long- slit spectra. The 
results are shown in Tables [2] and [3l 

Regions Rl and R2 are similar in physical conditions as well 
as in ionic abundances. In both regions, the abundances of the 
three species of oxygen as well as He^, Ar^^, N^, S^, and Cl^^ 
are enhanced with respect to the other regions. The opposite oc- 
curs for Ar^^ and He^^. 

With reference to elemental abundances, these are consistent 
with the mean values for PNe and Type I PNe, within uncer- 
tainties (see Table O. The exception is N, which show a defi- 
ciency, more in agreement with less-evolved objects, but con- 
sistent with previous estimates (Armour & Kingsburgh I2Q01I) . 
Given the values of these abundances, NGC 6309 would be a 
Type II PN (intermediate population), according to the classi- 
fication of Peimbert (1978 ). This is reinforced by the low N/0 
ratio, which is in agreement with the expected value for a low- 
mass central star (0.6 Mq, according to Kwok 2000). If this in- 
terpretation of the abundances is right, NGC 6309 would repre- 
sent another case of a non-Type I PN with bipolar outffows (see 
Vazquez et al. |1999al) . This strengthens the idea pointed out by 
Vazquez et al. that some low-mass central stars may also develop 
bipolar morphologies, in contrast to the fact that bipolar PNe are 
usually associated with relatively massive central stars (Corradi 
& Schwarz[T995l). 

2.5. Long slit high-dispersion spectroscopy 

We obtained high-dispersion optical spectra with the Manchester 
Echelle Spectrometer (MES) in the 2.1 m telescope (f/7.5) at the 
San Pedro Martir observatory (OAN-UNAM) during the period 
of 2001 May 22-23. We used a CCD SITe with 1024 x 1024 
pixels was used as detector, and set the slit width to 150-yum 
(1.6'0. A 2 X 2 binning was used, resulting in a spatial scale 
of Off 6 pixel'^ and a spectral scale of 0.1 Apixel'^ We centered 
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the spectral range at the Ha emission line. The spectra were 
wavelength calibrated with a Th-Ar arc lamp to an accuracy of 
±2kms"^. The achieved spectral resolution, as indicated by the 
FWHM of the comparison lines, was 12kms"^ Several spectra 
were obtained with the slit oriented E-W (slits A to F on Fig. [5]). 

We obtained on 2004 July 29-30, a second series of spectra, 
using the same telescope and instrument. In this case, the spectra 
were centered on the [O iii] /15007 emission line (slits G, H, I, 
and K on Fig. [5]), and the spectral scale was 0.08 Apixel'^ 

Position- velocity (PV) maps for the Ha and [O iii] A5001 are 
shown in Fig. [6l We note that the Heii/16560 line is also ob- 
served in the PV maps at the slits B, C, and D. In the PV maps 
(Fig.O relative position is measured with respect to the central 
star whereas radial velocity is measured relative to the systemic 
velocity, for which we deduce Vlsr = -32 ± 2kms"^ based on 
the expansion velocity at the position of the central star ([Oiii] 
spectra from slits H and K). This value is close to that reported 
by Schneider et al. (1983) of Vlsr ^ -33.4 ± 2.8 km s'^. 

In a first step to analyse the high-resolution spectra, we ex- 
tract from the PV maps the spatial position and radial velocity 
of the individual regions labeled in Fig.O The results are listed 
in Table (H The knots in the arms present a noticeable point- 
symmetry, both in their position with respect to the central star 
and in their radial velocity. In the case of the regions in the bright 
ellipse, Rl and R2 are moving at the systemic velocity while R3 
and R4, located along the minor axis, present the maximum ra- 
dial velocity in this structure. 

In a second step, we analyse the emission-line features in the 
PV maps as a whole. Slits H and K are tracing the minor and ma- 
jor axes of the internal ellipse, respectively. In the PV map of slit 
H (Fig. O, the regions E2 and W2 can be noted as bright con- 
densations located at the extremes of the emission-line features. 
Maximum expansion velocity at the lobes can be measured as 
2vexp - 80kms"^ (SW lobe expansion is best viewed in panel 
I). The lobe NE shows a bulk radial velocity of ^ -i-20kms"^ 
(as measured at 20'' from the center, panel H), whereas the lobe 
SW presents ^ -20kms"^ (as measured at 20'' from the center, 
panel I). The central region also seems to be expanding, showing 
a maximum velocity splitting of 2vexp - 69 ± 2kms"\ making 
it difl&cult to determine the radial velocity of the microstructures 
R3 and R4 accurately. However, the PV map from slit H ([O iii]) 
shows evidence of two additional components separated ^8" 
and ^ 46kms"^ that could be these features (locations are indi- 
cated on Fig. 6 by a pair of arrows). 

Slit K crosses the major axis of the ellipse. In the correspond- 
ing PV map, a velocity ellipse is observed. Features Rl and R2 
move with the systemic velocity. However, small regions beyond 
the ellipse are also found with radial velocities of ^ ±30 km s"^ . 
Maximum expansion velocity is observed at the center of the el- 
lipse (2vexp - 66 ± 2kms"^), which is nearly compatible with 
previous spectroscopic studies (Sabbadin I 19841 ^68 + 3 km s"^ ; 
Armour & Kingsburgh'SOOl , ^ 64 ± 5 km s"^). 

From these slits (H and K), the presence of two expanding 
asymmetrical lobes and a velocity ellipse in the center is evident. 
The expansion of these bipolar lobes is confirmed with data from 
the other slits, and the kinematics of their borders is traced by 
that of the point- symmetric knots. 

In the PV map from slit A, we detect an intense emission 
centered at 2" toward the west and relative velocity centered at 
Vr - -h20kms"^ This feature corresponds to the upper border 
of the ellipse (^ 5" above Rl). Faint emission from the halo 
is detected to the west of this feature (up to ^ 10" - 25") in 
agreement with the observed size of this structure. 



Table 4. Basic kinematical data for the main morphological fea- 
tures in NGC 6309. Angular distance and position angle (PA) are 
measured from the central star. Radial velocity is measured with 
respect to the systemic velocity (-32kms"^) from the emission 
line indicated in the last column. Uncertainties for the first row 
apply to the rest. 



Feature 


PA 

n 


Angular Radial 
distance velocity 
n (kms-i) 


Emission 
line 


Rl 


-15+4° 


9 ±r 


0+2 km s"^ 


[Oiii] 


R2 


+ 165 


10 





[Oiii] 


R3 


+75 


4 


+23 


[Oiii] 


R4 


-105 


4 


-23 


[Oiii] 


El 


+50 


30 


+35 


[Oiii] 


E2 


+60 


26 


+36 


Ha, [Oiii] 


E3 


+70 


27 


+35 


Ha 


E4 


+90 


25 


+24 


Ha 


Wl 


-130 


27 


-34 


Ha 


W2 


-120 


25 


-36 


Ha, [Oiii] 


W3 


-110 


24 


-34 


Ha, [Oiii] 


W4 


-90 


21 


-20 


Ha 



On the other hand, the PV map of slit I shows the radial 
velocity of feature E4 (the brightest feature) as ^ -h24kms"^ 
in addition to the expansion velocity of lobe SW close to ^ 
40kms"^ Finally, the PV map of slit C covers the region at 
lobe NE corresponding to the cone-like structure. This region 
expands to ^ 25kms"^ and is located from 5" to 15" toward 
the East from the center. The kinematics of the cone is practi- 
cally indistinguishable from that of lobe NE. 

3. Discussion 

3.1. The physical structure of NGC 6309 

The most reasonable explanation for the morphology and kine- 
matics of the central ellipse is that such structure corresponds to 
a tilted expanding torus. Although other possible structures are 
possible, this interpretation appears to be a better fit to our data. 
Assuming that the torus is circular and its diameter is ^ 20", 
we deduce an inclination angle of the torus axis with respect to 
the line of sight of 66°. The expansion velocity of the torus, cor- 
rected by the inclination angle is Vexp - 25kms"^. Therefore, 
the size of this structure and its kinematic age are ^ 0.2 pc 
and 3800 yr respectively, assuming a distance of 2 kpc. In addi- 
tion, the northeastern half of the torus is blueshifted whereas the 
southwestern half is redshifted. We note that these calculations 
have been obtained from the earth-based observations. Since the 
HST image shows a more complex structure in the torus, these 
results should be considered as an approximation only. 

With regard to the bipolar lobes, it is clear that a model 
of a single bipolar system with a main axis cannot reproduce 
the observed morphology, given the apparent two directions ob- 
served in the nebula. Therefore, we have considered a quadrupo- 
lar model consisting of two pairs of bipolar lobes, one of them 
oriented at PA -h76°, coincident with the torus axis, and the other 
at PA -h40°, coinciding with the apparent protrusions in the lobes 
along this axis. Each of these bipolar systems is described fol- 
lowing the formulation by Solf & Ulrich (.1985j , 

y(0) = Ve + (Vp - Ve) X sin"(|0|), (1) 
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Fig. 6. Position- velocity gray/contour maps corresponding to the slits A to F in the light of Ha (left panels), and slits G, H, I, and 
K in [O in] iSOO? (right panels). Slit orientation is indicated, being the same for all the left panels. In the right panels orientation 
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two features explained in the text. 
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where is the latitude angle above the equator; Ve and Vp are the 
equatorial and polar expansion velocity respectively; and a is an 
exponent that fits the specific hour-glass shape. With the main as- 
sumption that the axis of the torus and that of the bipolar ejection 
at PA -h76° are the same, the kinematics and morphology of this 
outflow can be reasonably reproduced with Ve = 25 km s"^ (ex- 
pansion velocity of the torus), Vp = 75 km s"^ and a = 6. There 
is no clue about the inclination angle of the outflow at PA -\-40° . 
If we assume the same inclination angle, we find reasonable fits 
for this bipolar system with Ve = 29kms"^ Vp = 86kms"^ 
and a = 6. In Fig. [71 we show the hour-glass shapes deduced 
from these fits overimposed on the unsharp-masking image of 
NGC 6309. The observed shape, as well as the kinematics of 
the knots, are reasonably reproduced by the models. With the 
model data and the assumption of a distance of 2 kpc, we find 
a kinematical age of 4 000 yr for the bipolar lobes at PA -h40° 
and of 3700 yr for the bipolar system at PA-h76°. These ages 
are not very diff'erent from each other and are also similar to the 
age deduced for the torus. This suggests that the formation of 
the diff'erent structures in NGC 6309 has occurred in a relatively 
small time span. In addition, as these kind of outflows have been 
observed in proto-planetary nebulae (PPN, e.g. Sahai 1998 ), our 
results impose a minimal value of time 4000 yr) for the PPN 
to PN stage in this nebula. Other structures such as the "cone" 
and the halo, are present in this nebula, but we must conduct 
further studies to understand their physical properties. 

3.2. The formation of NGC 6309 

The appearance of NGC 6309, as well as the position of the 
point- symmetric knots perpendicular to a radial position vec- 
tor from the central star, strongly resemble the point- symmetric 
structures observed in Fleming 1 (Lopez et al. 1993). Clifl'e et 
al. (1995 ) studied the interaction of a precessing jet with the 
interstellar medium. They find that point- symmetric structures 
can be formed by this interaction, and with the time, the indi- 
vidual bow- shocks can merge into a single shock structure. The 
ulterior evolution of such structures will lead to the formation 
of bipolar lobes with non-uniform brightness distribution being 
point- symmetric respect to the central star. On the other hand, 
Garcia-Segura & Lopez ( 2000 ) made models of bipolar PN with 
point- symmetric structures using a steady misalignment of the 
magnetic collimation axis with respect to the symmetry axis of 
the bipolar outflow. In their models, such morphology is also 
produced by the action of jets. 

A possible scenario for the formation of NGC 6309 is that 
the knots are the remnants of high- velocity bipolar collimated 
outflows that were ejected in the proto-planetary nebula (PPN) 
stage. There is evidence that precessing jets in young proto- 
PN carve bipolar or multipolar cavities in the initial st age of 
PNe evolution (e.g. Sahai & Trauger 1998^, Sahai et al. [20051 
Sanchez-Contreras et al. 2006). In such a scenario, the precess- 
ing jet is ejected during the proto-PN phase and after that, subse- 
quently, it is stopped and cooled by its surroundings. The knots 
follow the thermal expansion of the bipolar lobes. In this case, 
the brightness of the knots comes mainly from the photoionisa- 
tion process, as the shocked-cooled emission has turned ofl' or it 
is hidden by the efl'ect of a high-excitation photon source. This 
could explain the lack of [N ii] emission in the point- symmetric 
knots. 

In addition, we constructed a [Oiii]/Ha ratio image from 
images in Fig. [H Such ratio image can be used as a diagnos- 
tic tool, as has been shown by Medina et al. (120071 . Figure [8] 
shows the [O iii]/Ha ratio image in which an enhancement of the 



Fig. 7. NGC 6309 and the two bipolar outflows model. The fit 
of two hour-glass models is overimposed on the [O iii] unsharp- 
masking image presented before. The central torus is also drawn. 




Fig. 8. [O iii]/HQf ratio image. White represents high values of 
the ratio. 



[O III] emission is seen in the edge of the lobes, clearly related 
to the point- symmetric knots. It is noticeable that the emission 
vanished at the end of the SW lobe, possibly corresponding to a 
breaking of the s hell du e to the action of the outflow. According 
to Medina et al. (120071) . the [Oiii] enhancement would be com- 
patible with structures shaped by collimated outflows, as in the 
case of IC 4634. Finally, the circular halo probably corresponds 
to the remnant of the envelope ejected as a slow wind from the 
central star when it was at its AGB phase. 
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4. Conclusions 

We have carried out an analysis of the planetary nebulae 
NGC 6309, based on ground-based and space-based imaging, 
high- and low-dispersion spectroscopy, as well as VLA-D radio 
continuum. We summarize the main conclusions of this work as 
follows. 

NGC 6309 can be described as a quadrupolar PN formed by a 
bright central torus, two systems of bipolar lobes oriented at dif- 
ferent directions, and point- symmetric knots that trace the edges 
of the lobes. The torus expands at 25kms"^ whereas the po- 
lar expansion velocity of the lobes is 75 km s"^ for a first bipolar 
system at PA 76°, and 86 km s"^ for another bipolar system at PA 
40°. Assuming a distance of 2kpc for the nebula, the kinematic 
ages of the structures ranges from 3700 to 4000 yr, suggesting 
that they have been formed in a short time span. The knots at the 
edges of the lobes suggest that the lobes have been formed by 
rapidly precessing bipolar jets that carved cavities in the previ- 
ous red giant envelope. In addition to these structures, we detect 
a circular halo surrounding the torus and bipolar lobes, which 
probably corresponds to the envelope ejected in the AGB phase 
by the central star. There is also a conelike structure embedded 
in one of the lobes and with its base on the torus. 

We also study internal variations of the physical conditions 
and chemical abundances in NGC 6309. The low-dispersion 
spectra indicate a high-excitation nebula, with low to medium 
variations of its internal physical conditions (10, 600 K < 
re[Oiii] < 10, 900 K; 10, lOOK < T^Nii] < 11, 800 K; 
1440cm-^ < A^e[Sii] < 4000cm-^; HOOcm"^ < A^e[Cliii]< 
2600 cm"^; 1000 cm"^ < A^e[Ariv] < 1700 cm"^). The ra- 
dio continuum emission indicates a mean electron density of 
^ 1900cm"^; emission measure of 5.1 x lO^pccm"^; and an 
ionised mass M(Hii)^ 0.07 M©. The logarithmic extinction co- 
efficient cufi range from 0.70 to 0.97 (mean value 0.87), which 
probably is produced by differences in the internal dust distribu- 
tion. 
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Table 1. Dereddened spectra and physical conditions from different regions in NGC 6309. Intensity values are relative to /^s=100. 
Extinction law (f^) is from Howarth (1983). The size of the integration region (in arcsec) along the slit is 7! 4 for Rl, R2, R3, and 
R4, and 5'.'3 for the rest. Slit width was set at Vf6 in all the cases. 



Ion 


^(A) 




Rl 


R2 


R3 


R4 


E2 


Wl 


W2 


[On] 


3727 


0.256 


108 


102 












[Ne III] 


3869 


0.230 


163 


165 












[Ne III] 


3968 


0.210 


68 


67 












[S II] 


4071 


0.189 


6.5 














US + He II 


4101 


0.182 


33 


34 












Hy 


4340 


0.127 


55 


53 


26 


27 


52 


— 


— 


[Oiii] 


4363 


0.121 


12 


12 


7.9 


8.5 


— 


— 


— 


Hex 


4471 


0.095 


8.2 


7.1 


3.4 


3.2 








Hen 


4540 


0.078 


— 


— 


2.8 


3.3 


— 


— 


— 


Nin 


4640 


0.054 


4.3 


5.0 


5.2 


6.1 








Hen 


4686 


0.043 


29 


31 


81 


79 


86 


106 


— 


[ Ar iv] 


4711 


0.037 


6.0 


5.8 


12 


12 


— 


— 


— 


[ Ar iv] 


4740 


0.030 


4.6 


4.6 


9.4 


9.3 








H/3 


4861 


0.000 


100 


100 


100 


100 


100 


100 


100 


[Oiii] 


4959 


-0.024 


487 


464 


335 


365 


341 


364 


451 


[Oiii] 


5007 


-0.036 


1504 


1466 


993 


1081 


1046 


1097 


1238 


Ni 


5199 


-0.082 


0.9 














Hen 


5411 


-0.133 


2.5 


2.3 


6.2 


6.0 


8.5 


— 


— 


[Clni] 


5517 


-0.154 


1.1 


1.0 


0.7 


0.7 


— 


— 


— 


[CI in] 


5537 


-0.157 


1.0 


1.0 


0.7 


0.6 


— 


— 


— 


[Nil] 


5755 


-0.195 


1.0 


1.0 


0.6 


— 


— 


— 


— 


He I 


5876 


-0.215 


13 


12 


7.1 


8.1 


— 


— 


— 


[Oi] 


6300 


-0.282 


6.8'' 




2.5 


1.2 


— 


— 


— 


[S m] 


6312 


-0.283 


— 


— 


1.8 


1.7 


— 


— 


— 


[Oi] 


6364 


-0.291 


1.9 


1.8 


0.8 


0.4 


— 


— 


— 


[Arv] 


6435 


-0.302 






1.4 


1.4 








[Nil] 


6548 


-0.318 


24 


31 


9.3 


4.7 








Ha 


6563 


-0.320 


285 


285 


283 


284 


284 


284 


284 


[Nn] 


6583 


-0.323 


60 


60 


27 


12 


8.8 






He I 


6678 


-0.336 


3.4 


3.6 


2.4 


2.6 


— 


— 


— 


[S II] 


D/i / 


-U.J4z 


o.U 






l.J 








[S III 


6731 


-0.344 


11 


12 


4.8 


2.1 








TArvl 

Li-vi vj 


7006 


-0.380 






3.1 


2.7 








He I 


7065 


-0.387 


3.1 


3.1 


1.9 


2.0 








[ Ar in] 


7136 


-0.396 


17 


19 


14 


13 


19 


21 


21 


[Any] + [Fen] 


7170 


-0.401 


1 




0.9 


0.3 








Hen 


7178 


-0.401 




— 


0.9 


0.7 


— 


— 


— 


[0 II] 


7320 


-0.419 


5.6° 


5.7^^ 


2.3 


1.4 








[On] 


7330 


-0.420 






2.0 


1.3 








CllJ3 






0.96 


0.97 


0.70 


0.80 


0.86 


0.85 


0.97 


log lufi 






-13.4 


-13.5 


-13.2 


-13.4 


-15.2 


-15.3 


-15.9 


T,[Om\ 


(K) 




10 800 


10900 


10700 


10600 








re[Nii] 


(K) 




10 300 


10100 


11800 










A^e[Sll] 


(cm-3) 




1800 


1440 


4000 


3600 








A^e[CllIl] 


(cm-3) 




2100 


2350 


2600 


1700 








A^e[Ariv] 


(cm-3) 




> 1000 


> 1500 


> 1700 


> 1500 








logHQf/[Nii] 






0.533 


0.499 


0.889 


1.244 








logHQ'/[S n] 






1.182 


1.142 


1.565 


1.913 








log[Oiii]5007/HQ' 






0.723 


0.711 


0.545 


0.580 








log^6717M6731 






0.748 


0.792 


0.603 


0.623 









^ This emission line is blended with [S iii]6312A. 
^ This emission line is blended with [O ii]7330A. 
Mean uncertainties are Acujs = ±0.02, AT^ = +500 K, AA^e 



= +400 cm-^ 
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Table 2. Ionic abundances of NGC 6309. In the case of ions with more than one transition, a flux- weighted average was performed. 
Electron temperature and density values of Tg =10700 K and A^^ = 2400 cm"^ were used on the calculations. 



Ion 






Rl 


R2 


R3 


R4 


E2^ 




W2^ 


He+2 






0.02 


0.03 


0.07 


0.06 


0.06 


0.09 




He+ 






0.09 


0.09 


0.05 


0.06 






— 




(xlQ- 




4.68 


4.93 


2.02 


3.78 


3.66 


3.82 


4.41 


0+ 


(xlQ- 




4.70 


5.93 


1.28 


2.16 










(xlQ- 




1.13 


1.29 


0.27 


0.25 


— 


— 


— 


Ar+4 


(xlQ- 




— 


— 


4.26 


5.83 


— 


— 


— 


Ar+3 


(xlQ- 




8.78 


9.22 


12.2 


19.4 








Ar+2 


(xlQ- 




1.48 


1.68 


0.90 


1.18 


1.71 


1.91 


1.92 


N+ 


(xlQ- 




1.15 


1.37 


0.37 


0.25 


0.18 






NO 


(xlQ- 




2.18 
















(xlQ- 








2.14 


3.99 










(xlQ- 




5.29 


5.87 


2.07 


1.31 








Ne+2 


(xlQ- 




1.59 


1.75 












Cl+2 


(xlQ- 




1.28 


1.28 


0.63 


0.92 









Table 3. Elemental abundances of NGC 6309. ICFs were obtained following Kingsburgh & Barlow (1994). Comparison with mean 
values of All-type PNe, Type I PNe and other objects are shown. Except for He, all the abundances relative to H are logarithmic 
values with H=-i-12. 



Ion 


Rl 


R2 


R3 


R4 


G04" 


PNe^ 


TIPNe'^ 


Hii^ 


Sun^ 


He/H 


0.11 +0.01 


0.11+0.01 


0.12 + 0.01 


0.12 + 0.01 


0.10 + 0.01 


0.12 + 0.02 


0.13 + 0.04 


0.10 + 0.01 


0.09 + 0.01 


0/H 


8.74 + 0.02 


8.72 + 0.02 


8.72 + 0.01 


8.74 + 0.01 


8.64 + 0.02 


8.68 + 0.15 


8.65 + 0.15 


8.70 + 0.04 


8.66 + 0.05 


N/H 


7.93 + 0.02 


7.92 + 0.02 


8.04 + 0.01 


7.89 + 0.02 


8.20 + 0.12 


8.35 + 0.25 


8.72 + 0.15 


7.57 + 0.04 


7.78 + 0.06 


Ne/H 


8.26 + 0.02 


8.25 + 0.02 






7.82 + 1.00 


8.09 + 0.15 


8.09 + 0.15 


7.90 + 0.10 


7.84 + 0.06 


Ar/H 


6.40 + 0.10 


6.42 + 0.10 


6.54 + 0.01 


6.51+0.01 




6.39 + 0.30 


6.42 + 0.30 


6.42 + 0.04 


6.18 + 0.08 


S/H 


6.74 + 0.01 


6.73 + 0.01 


6.84 + 0.01 


6.87 + 0.01 


6.49 + 0.05 


6.92 + 0.30 


6.91+0.30 


7.06 + 0.06 


7.14 + 0.05 



^ Taken from Gornv et al. ([20041) . 

^ Average for PNe (Kingsburgh & Barlow [T9941 ). 

Average for Type I PNe (Kingsburgh & Barlow [T9941 ). 
^ Average for H ii regions (Shaver et al. 198T). 
^ The Sun (Grevesse, Asplund, & Sauval 2007) . 



